Equation of state of cosmic strings with fermionic current-carriers 
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The relevant characteristic features, including energy per unit length and tension, of a cosmic 
string carrying massless fermionic currents in the framework of the Witten model in the neutral limit 
are derived through quantization of the spinor fields along the string. The construction of a Fock 
space is performed by means of a separation between longitudinal modes and the so-called transverse 
zero energy solutions of the Dirac equation in the vortex. As a result, quantization leads to a set 
of naturally defined state parameters which are the number densities of particles and anti-particles 
trapped in the cosmic string. It is seen that the usual one-parameter formalism for describing the 
macroscopic dynamics of current-carrying vortices is not sufficient in the case of fermionic carriers. 

98.80.Cq, 11.27.+d 



I. INTRODUCTION 



The mechanism of spontaneous symmetry breaking involved in early universe phase transitions in some Grand 
Unified Theories (GUT) might lead to the formation of topological defects [jlj. Among them, only cosmic strings 
happen to be compatible with observational cosmology if they form at the GUT scale. It was shown however by 
Witten [0] that, depending on the explicit realization of the symmetry breaking scheme as well as on the various 
particle couplings, a current could build along the strings, thereby effectively turning them into superconducting 
' wires. Such wires were originally considered in the case the current couples to the electromagnetic field so they may 
be responsible for a variety of new effects, including an explosive scenario for large scale structure formation for 
CO 1 which an enormous energy release was realized in the form of an expanding shell of non propagating photons in the 
surrounding plasma || . 

The cosmology of strings has been the subject of intense work in the last twenty years or so Q, mainly based 
O |' on ordinary strings, global or local, aiming at deriving the large scale structure properties stemming from their 
JL , distribution as well as their imprint in the microwave background It was even shown Q that the most recent 
data |Q] might support a non negligible contribution of such defects. As such a result requires ordinary strings, it 
turns out to be of uttermost importance to understand the influence of currents in the cosmological context. 

Indeed, it can be argued that currents might drastically modify the cosmological evolution of a string network: The 
most clearly defined consequence of the existence of a current flowing along a vortex is the breaking of the boost 
invariance, since the current itself defines a privileged frame. In other words, the energy per unit length U and the 
tension T become two different numbers, contrary to the ordinary (Goto-Nambu |^]) case. As a result, string loops 
become endowed with the capability of rotation (the latter being meaningless for U = T), and the induced centrifugal 
force permits equilibrium configurations, called vortons ||. They would very rapidly reach a regime where they would 
scale as ordinary non relativistic matter, until they come to completely dominate the Universe fl(i(| . 

In the original Witten model currents could form by means of two different mechanisms. Scalar fields, directly 
coupled with the string forming Higgs field, could feel a localized potential into which they could accumulate in 
the form of bound states, while fermions could be trapped along the string, propagating at the speed of light, as 
zero energy solutions of the two dimensional Dirac equation around the vortex. Other models were proposed where 
fermions could also propagate in the string core at lower velocities in the form of massive modes [jll| , or (possibly 
charged) vector fields could also condense All these models have essentially made clear that the existence of 

currents in string is much more than a mere possibility but rather an almost unavoidable fact in realistic particle 
physics theories. 

For scalar as well as vector carriers, the task of understanding the microphysics is made simple thanks to their 
bosonic nature: all the trapped particles go into the same lowest accessible energy state and the field can be treated 
classically [Jl3|,[l4j . Even the surrounding electromagnetic jl5| and gravitational [pj6|p^7| fields generated by the current 
can be treated this way and the backreaction can be included easily [l§| ]. 

Meanwhile, a general formalism was set up by Carter [l9| to describe current-carrying string dynamics. The 
formalism is based on a single so-called state parameter, w say, of which the energy per unit length, the tension 
and the current itself are functions. Such a formalism relied heavily on the fact that for a bosonic carrier, the 
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relevant quantity whose variation along the string leads to a current is its phase, and the state parameter is essentially 
identifiable to this phase gradient. Various equations of state relating the tension to the energy per unit length were 
then derived (^0|, based on numerical results and the existence of a phase frequency threshold [O. It even includes 
the special case of a chiral current [EM, although the latter originates in principle only for a purely fermionic current. 
Therefore, it was until now implicitly assumed that such a formalism would be sufficient to describe whatever current- 
carrying string configuration. It will be shown in this paper that this is in fact not true and an extended version, 
including more than one parameter, is needed p^ |. 

The state parameter formalism, apart from being irrelevant for fermionic current-carrying strings, can only provide 
a purely classical description of their dynamics. This is unfortunate since the most relevant prediction of supercon- 
ducting cosmic string models in cosmology is the existence of the vorton states discussed above. These equilibrium 
configurations of rotating loops are not necessarily stable, and in fact, this is perhaps the most important question 
to be answered on this topic. Indeed, any theory leading to the the existence of absolutely stable vortons predicts 
a cosmological catastrophe jlO| and must be ruled out. One may therefore end up with a very stringent constraint 
on particle physics extension of the standard model of electroweak and strong interactions. To decide clearly on this 
point requires to investigate both the classical and the quantum stability of vortons. 

Classical stability has already been established in the case of bosonic carriers [^3| for whatever equation of state |24|] 
on the basis of the one parameter formalism. Yet it will also have to be addressed in the more general context that will 
be discussed below. In the meantime it was believed that a quantum treatment was necessary in order to decide on 
the quantum stability: as one wants to compare the characteristic life-time of a vorton with the age of the Universe, 
quantum effects can turn out to be relevant; hence the following work in which the simplest of all fermionic Witten 
models is detailed that can give rise to both spacelike as well as timelike charge currents, generalizing the usual 
point of view p5|] . 

Let us sketch the lines along which this work is made. 

A two-dimensional quantization of the spinor fields involved along a string is performed. Owing to anti-particle 
exitation states, one can derive the conditions under which the current is of arbitrary kind. Moreover, an equation of 
state giving the energy per unit length and the tension is obtained that involves four different state parameters which 
are found to be the number densities of fermions, although three of them only happen to be independent. 

In section ||, the model is presented and motivated, and the equations of motion are derived. Then in section III , 
we obtain plane wave solutions along the string by separating transverse and longitudinal dependencies of spinor fields 
in the vortex. The zero mode transverse solutions are then constraints to be normalizable in order to represent well 



defined wave functions. The quantization restricted to massless longitudinal modes is performed in section IV. As a 



result, the classical conserved currents obtained from Noether theorem, like energy- momentum tensor and fermionic 
currents, are expressed in their quantum form. All these quantum operators end up being functions of the fermionic 
occupation numbers only. In the last section (section [v|) , the classical expressions for the energy per unit length and 
the tension are derived and discussed from computation of quantum observable values of the stress tensor operator 
in the classical limit. Contrary to the bosonic current-carrier case where there is only one state parameter jlj], the 
classical limit of the model involves four state parameters in order to fully determine the energy per unit length and 
the tension. The cosmological consequences of this new analysis are briefly discussed in the concluding section. 



II. EQUATIONS OF MOTION 



We are going to be interested in the purely dynamical effects a fermionic current flowing along a cosmic string may 
have. The model we will be dealing with here is a simplified version of that proposed by Witten which involves two 
kinds of fermions, in the neutral limit. This limit, for which the coupling between fermions and electromagnetic-like 
external fields is made to vanish, permits an easy recognition of the dynamical effects of the existence of an internal 
structure as in Ref. [Q. 



A. Particle content 



The model we shall consider involves a complex scalar Higgs field, <E> say, with conserved charge qc$ under a local 
t/(l) symmetry, together with the associated gauge vector field B^. In this simple Abelian Higgs model vortices 
can form after spontaneous breaking of the U(l) symmetry. The minimal anomaly free model |^] with spinor fields 
requires two Dirac fermions denoted \t and X ', with opposite electromagnetic-like charges, getting their masses from 
chiral coupling with the Higgs field and its complex conjugate. They also have conserved gauge charges from invariance 
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under the broken U(l) symmetry, qc^ R1 qc^ L , qc XR , and qc XL , for the right- and left-handed parts of the two fermions 
respectively. The Lagrangian of the model therefore reads 

C = C h + C g + + L x (1) 

with Lh, C g and C^, C x , respectively the Lagrangian in the Higgs, gauge, and fermionic sectors. In terms of the 
underlying fields, they are 

£ & = (£>"*) (2) 



C g = --H^H^, (3) 



C x = l - \X-fD»X - (D^X^X] - gX±±^X<t>* - gX^—J^X^, (5) 
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where we have used the notation 

(6) 

Br+iq **- ** ^)*, (7) 
B.+iq ^'^ ^B^X, (8) 

(9) 
(10) 

The equivalence with the Witten model |2j appears through a separation into left- and right-handed spinors. Let us 
define ^>r and ^l, respectively the right- and left-handed parts of the Dirac spinor field VP (and the same for X), 
eigenvectors of 75, 

*«=^-y^*, and * L = i^p*. (11) 
The Lagrangian for the spinor field \I/ now reads 

A/, = \ ^rI^D^r - (D^ R )^ R ] + l - ^ L y i D^ L - (D^^^l] - g¥ L * R $ - 5 W L $*, (12) 
with the associated covariant derivatives 

D^r(l) = (V M + iqc^ R(L) B^ R{L) . (13) 

It is clear with the Lagrangian expressed in this way that the invariance of the action under U(l) transformations 
requires 

C 1pL ~ C 1pR = C 4> = C XR ~ c xl- (14) 

B. Equations of motion 

As we wish to deal with a cosmic string, the Higgs and gauge fields can be set as a vortex-like Nielsen-Olesen 
solution and they can be written in cylindrical coordinates as follows pTj| 

$ = ^(r)e fa W, B tt =B(r)5 fl8 . (15) 

In order for the Higgs field to be well defined by rotation around the string, its phase has to be proportional to the 
orthoradial coordinate, a(6) = n6, where the integer n is the winding number. The new fields tp and a are now real 
scalar fields and are solutions of the equations of motion 
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where 



Qju = V M a + qc^Bft, 

m^p = gip cos a + igipjs sin a, 

m x = cos a — ig<f"f5 sin a. 

In the same way, the equations of motion for the gauge and spinor fields are 



i (V M tf) 7 M 



97 M _ 



dX 



(16) 
(17) 



(18) 
(19) 
(20) 



(21) 
(22) 

(23) 
(24) 
(25) 



The fermionic currents j^( x \ have axial and vectorial components due to the different coupling between left- and 
right-handed spinors to the gauge field. The two kinds of current are required to respect gauge invariance of the 
Lagrangian. In terms of spinor fields, they read 



(26) 



with 



•LI 

3 ipv 



'Xv 



= ^ 



-XR ~ ^Xl 

2 

-Xr ~ c Xl ■ 



(27) 



III. TRANSVERSE SOLUTIONS AS ZERO MODES 



A. Plane wave solutions 



The study of the fermionic fields trapped along the string can be performed by separating the longitudinal and 
transverse solutions of the equations of motion. The plane wave solutions are therefore expressed in the generic form 



4 ±: 



_±£(o>f— kz) 



/ £i(r)e-"™ ie \ 
\ \ 4 {r)e- imie J 



X, 



(±) _ „±i(urt-fez) 



/Ci(r)e- ai9 \ 
( 3 (r)e-^ e 



V C 4 (r)e 



-u 4 e 



(28) 



Inside the vortex, the numbers m, and k have to be integers in order to produce well defined spinors by rotation 
around the string. In the following, the Dirac spinors will be expressed in the chiral representation, and the metric is 
assumed to have the signature (+,—,—,—). Plugging the expression (|8|) into the equations of motion ( p2|) and j2^ ) 
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yields the differential system 

+ ~ [-<l c ^R B ( r ) + m i] ftO")) e 



-i(mi-l)6 



igip(r)e 



— i(m.4-\-n)0 



dr r 



3 -i(m 2 + l)0 _ 



'■■g<fi(r)e 



^ + - l-qc^ L B(r) + ms] &(r) ) e^™*" 1 )* + w^e"^""^^) 



T«(fc- w)C4(r)e-"" 4e , 
±i(fc + w)C 3 (r)e- im3e . 



(29) 



and 



Similar equations are obtained for the field X with the following transformations, £ — * £, c^, R{L) — > c XH(t) , 
n — > — n. because of its coupling to the anti-vortex instead of the vortex. Note that if, instead of the vectorial phases 
ansatz (pq), we had chosen a matricial phases ansatz in the form 



\J/(±) — e ±i(wt-fez) 



/ £n(r)e- imilfl + . . . + CwWe"™ 149 \ 
&i(r)e~' im2ie + • • ■ + £ 2 4 Me"™ 249 
6i(r)e" mi319 + . . . + C 3 4(r)e- im3 ^ 

\ ^4i(r)e- im4lff + . . . + e 4 4(r)e- jm44e / 



(30) 



we would have found, from the requirement of having at most four independent phases in the equations of motion, 
that the matrix m,j has to verify for all i, rriij = rriik, for all (j,k). Consequently, the vectorial ansatz ( p8| ) is the 
most general for solutions with separated variables. 



B. Transverse solutions 



From the differential system ( fj9|) , it is obvious that the four phases m, cannot be independent parameters if the 
spinors fields are not identically zero. It is also impossible to find three independent phase parameters since each 
equation involves precisely three different angular dependencies. The only allowed angular separation requires two 
degrees of freedom in 0, and the only relevant relation for trapped modes in the string reads from Eq. ( p9[ ) 

mi — 1 = 7714 + n, rri2 + 1 = + n. (31) 



1. Zero modes 



Introducing the two integer parameters p = mi and q = 7713, and using equation (|3l]), for p q + n, the system 
(pit) reduces to the set 



^ + - [qc ilR B(r) -(q + n- 1)] &(r) - iff¥>(r)&(r) = 0, 
dr r 

^ + - hqc^ L B(r) + q] &(r) + ig<p(r)&(r) = 0, 
dr r 

(fc + w)6(r) = 0, 
I (fc + w)6(r) =0, 



(32) 



-jT + - [-?c^ R B(r) +p] £i(r) - igip(r)^(r) = 0, 

^ + i [gc^S(r) - (p - n- 1)] &(r) + w(r)6M = 0, 

(k-w)£i(r) = 0, 
I (fc-w)£4(r) = 0. 



(33) 



There are two kinds of solutions which propagate along the two directions of the string at the speed of light and which 
were originally found by Witten [g]: Either fc = uj and £2 = £3 = 0, or k — —lu and £1 = £4 = 0. These zero modes 
must also be normalizable in the transverse plane of the string in order to be acceptable as wave functions. 
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2. Index theorem 



In the two cases k — to and k = —uj, we will call the corresponding zero modes, X^,(r, 9) and Y^(r, 
of the systems (H|) and (|33"1), respectively, i.e., 



/o 



Y h 



re 



/a 






-i(p- 



the solutions 



(34) 



These have to be normalizablc in the sense that J |X^,| 2 r dr dd and J {Y^r dr d9 must be finite. Thanks to the 
regularity of the vortex background, the divergences in these integrals can only arise close to the string core or 
asymptotically far away from it. As a result, it is sufficient to study asymptotic behaviors of the solutions to decide 
on their normalizability p7| . 

Let us focus on the zero mode X^ solution of Eq. ( |32| ) , keeping in mind that Y^ can be dealt with in the same way. 
The asymptotic behaviors at infinity are easily found as solutions of the limit at infinity of the differential system 
( |32| ) . Note that the identification between the equivalent solutions and the solutions of the equivalent system is only 
allowed by the absence of singular point for the system at infinity, and it will not be so near the string since r — is a 
singular point, so that the Cauchy theorem does no longer apply. From Eq. (p^), the eigensolutions of the equivalent 
system at infinity are in the form exp (ztgryr), and thus, there is only one normalizable solution at infinity. 

Near the string, i.e., where r — > 0, the system is no longer well defined, the origin being a singular point. Approximate 
solutions can however be found by looking at the leading term of a power-law expansion of both system and functions, 
as originally suggested by Jackiw and Rossi Because the Cauchy theorem does no longer apply, many singular 
solutions might be found at the origin, and among them, the two generic ones which match with the two exponentials 
at infinity. Near the origin, the Higgs and gauge fields are known to behave like J3(J 



ip(r) ~ ip r 



\n\ 



B(r) oc r 2 



so that the leading contribution near the string of the zero modes can be found as 
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6(r) 



a 2 r 
a 3 r° 



(35) 



(36) 
(37) 



with Oi and on real parameters to be determined. The values of the exponents on are therefore given by the leading 
order terms in system (p2) , and one obtains three solutions, the first one of which being singular, 



a 2 r 



q+n— 1 



provided 



< q < 1 



a 3 r-i ) 2*2 
The two other solutions are the generic ones which have to match with the solutions at infinity 
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gl 



-q+M + l 



a 3 r" 
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q+n— 1 



a 3 r' 



q+\n\+n 



(38) 



(39) 



where the relationships between the parameters at have not been written, since they are clearly obtained from Eq. ( J3S 
Normalizability near the origin requires that the integrals J |^2| 2 rdr and / j^j 2 rdr converge, and this yields 



— n < q < 1. 

Analogous considerations for the system ( ^3|) show the convergence criterion in this case to be 

n < p < 1. 



(40) 



(41) 



The number of sets of parameters p and q satisfying the previous inequalities is precisely the number of well defined 
zero modes, respectively Y^ and X^, which are also normalizable. In order to match with the single well behaved 
solution at infinity, the two independent solutions near the string have to be integrable. Therefore, for a vortex solution 
with a positive winding number n, there are only n normalizable zero modes, which are the X^ ones. Similarly in 
the case of an anti- vortex with negative winding number — |n|, one finds also |n| zero modes Y^. This is the index 



theorem found by Jackiw and Rossi 



Recall that the model involves two kinds of fcrmions, and all the previous 
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considerations apply as well for the field X with the simple transformation n — > — n. Therefore, the normalizable zero 
modes are swapped compared to those of the field 

Finally, for a vortex with positive winding number n, there are always n massless plane wave solutions for both 
spinor fields, which read 



p 



3 ±ifc(t+z) 



±ifc(i-z) 



/o 

^ 2 ( r ) e -«(9+™-l)S 

C 3 (r)e-^ 
V° 

/Ci(r)e"^ 



V C4(r)e- i ( p+n - 1 ) 9 



5±l fc(t+-)^( r , 0), 



with now q = m 3 and p = h which satisfy 

—n < q < 1, and — n < p < 1. 
Note that they are eigenvectors of the j°j 3 operator, and they basically verify 
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= 0, 
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-16 
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The X zero modes X x and Y x verify the same relationships with replaced by £i 



5. Massive modes 



(42) 



(43) 



(44) 



(45) 



The case mi — m 3 + n allows four-dimensional solutions of the system (|29|). In Particular, these solutions do no 
longer require u> — ±fc and therefore represent massive modes. As before, the interesting behaviors of these modes 
are found by studying the solutions of the equivalent system asymptotically and by looking for the leading term of a 
power-law expansion of both system and solution near the string core. 

At infinity, the system is well defined and there are two twice degenerate eigensolutions exp (±ilr) out of which two 
are normalizable, with 



SI = \/ g 2 f] 2 — (uj 2 — k 2 ). 



(46) 



Near the origin, at r = 0, the system is singular, and because of its four dimensions there are much more singular 
solutions than the previous two dimensional case, and among them the four generic ones which match with the four 
ones at infinity. The leading term in asymptotic expansion can be written in a standard way 



Si{r)~a<r*. (47) 

Plugging these expressions in the system ( p9| ) with Eq. (|3^), and keeping only leading terms at r = gives, after some 
algebra, the four generic solutions 
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with m = mi, and where the relationships between the coefficients have not been written as they are essentially 



given by a linear system in at given by Eq. (|29|). The solutions (48) will be normalizable near the string if, for all 
i, J |£i| 2 rdr is hnite. Moreover there will be, at least, always one massive bound state if there are at least three 
normalizable eigensolutions to match with the well-defined ones at infinity. This is only allowed if the parameter m 
verifies simultaneously three of the following conditions 

sup(0,n) <m<mf (1,1 + n). (49) 

Because m is necessary an integer, this condition cannot be achieved. This criterion, originally derived and used by 
Jackiw and Rossi in order to enumerate the number of zero modes in a vortex- fermion system ]27| , is only sufficient 
and thus, normalizable massive bound states may exist, but are model dependent since it is necessary that a particular 
combination of the two normalizable eigenmodes near the string core match exactly with a particular combination of 
the two well-defined ones at infinity. 

Such massive bound states depend therefore of the particular values of the model parameters. Recently, it was shown 
numerically [^5| that the Abelian Higgs model with one Weyl fermion admits always at least two massive bound states, 
as a result, the present toy model also may have such states. However, in order to simplify the quantization, we will 
only consider the generic zero modes, and consequently, the following results will be relevant for cosmic string only 
when the occupancy of the massive bound states can be neglected compared to the occupancy of the zero mode states. 
Such physical situations are likely to occur far below the energy scale where the string was formed, since the massive 
states are generally expected to decay much more rapidly than the massless ones f25| . 

The generic massless normalizable transverse solutions of the fermionic equations of motion in the string with 
winding number n are the n zero modes. For the spinor field VE" coupled with the vortex, we find that the particles 
and the anti-particles can only propagate at the speed of light in one direction, "— z" say direction along the string, 
whereas the spinor field X propagates in the opposite, direction. The existence of such plane waves allows us 

to quantize the spinor fields along the string. The zero modes themselves will therefore be transverse wave functions 
giving the probability density for finding a trapped mode at a chosen distance from the string core. 



IV. FOCK SPACE ALONG THE STRING 



The spinor fields can be expanded on the basis of the plane wave solutions computed above. A canonical quantization 
can then be performed along the z-axis which provides analytical expressions for these fields in two dimensions once 
the transverse degrees of freedom have been integrated over. It is therefore possible to compute the current operators 
as well as their observable values given by their averages in a particular Fock state. In the following we shall take 
a vortex with a unit winding number n — 1 and the subscript of the zero modes X^ and Y x will be forgotten since 
there is not possible confusion. 



A. Canonical quantization 



We shall first be looking for a physical expansion of the spinor fields in plane waves, in the sense that creation and 
annihilation operators are well defined. The Hamiltonian is then calculable as a function of these and will be required 
to be positive to yield a reasonable theory. 



1. Quantum fields 



As shown above, the spinor fields ^ and X propagate in only one direction, therefore in expressions (|42|) and (|42^ 
the momentum k can be chosen positive definite. Let us, once again, focus on the spinor field ^P. The natural way to 
expand it in plane waves of positive and negative energies is 



dfc 
2vr2fc 



&*(-*)■ 



s ife(t+z) 



b(-k)e 



-ik(t+z) 



X. 



(50) 



The Fourier transform of 'J on positive and negative energies has been written with similar notation b* and b unlike in 
the free spinor case. Indeed, note that, in the string, the zero modes are the same for both positive and negative energy 
waves, so that the only way to distinguish particles from anti-particles is in the sign of the energy. The integration 
measure dfc/(27r2fc) is the usual Lorentz invariant measure in two dimensions. Note that k is chosen always positive 
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in order to represent physical energy and momentum actually carried by the field along the string; hence the negative 
sign in b'(—k) and b(—k), which is a reminder that the spinor field W propagates in the "— z" direction. In the same 
way, the field X is expanded as 



X = 



/ 2?r2fc 

The Fourier transform will be written with the normalization convention 



Ah 

tft(k)e tt (*-") +d{k)e- ik ^ Y. (51) 



dze t(k- k ')z = 27T S(k- k'). (52) 



2. Creation and annihilation operators 

The Fourier coefficients can be expressed as functions of the spinor field \t or With equation ( |50| ) and $52\), let 
us compute the following integral 

r Ay \\X\\ 2 
rdrd9dze lk{ - t+ ^X^ = / — [b\-k')8{k + k') + b{-k')8(k - k')] ||X|| 2 = l ±-£-b(-k), (53) 

where we have defined 

\\X\\ 2 = f rdrd9\X\ 2 . (54) 



Note that the separation between b and b^ only arises from the chirality of the spinor field because the integration is 
performed only over positive values of the momentum k; this is why the S(k + k') term vanishes. In the following we 
will assume that the zero modes are normalized to unity, ||X|| 2 = 1, and \\Y\\ 2 = 1. Playing with similar integrals 
gives us the other expansion coefficients 



(55) 



(56) 



b(-k) = 2k J rdrd9dze lk( - t+z) X^ 7 ^(-k) = 2k J rdrdOdze- lk{t+z) ^X, 

b\-k) = 2k J rdrd6dze- tk{t+z) X^, b(-k) = 2k J r dr d8 dz e lk{t+z) ^ X , 
and the corresponding relations for the spinor field X 

*(») . »/,*d,d,.-<~>l**, m . »f rir »<,+0-4#r, 

d{k) = 2kj rdrdddze ik(t - z) Y^X, d\k) = 2k J r dr d6 dze^ 1 -^ X ] Y . 

From these, one gets the necessary relations to define creation and annihilation operators 

b\-k) = [fe(-fc)] f , and fct(-fc) = [b(-k)] 1 , 
dt(Jfe) = [d(k)] j , and d\k) = [d(fc)] f . 

3. Commutation relations 

The canonical quantization is performed by the transformation of Poisson brackets into anticommutators. Here, we 
want to quantize the spinor fields only along the string, and therefore let us postulate the anticommutation rules at 
equal times for the quantum fields 

{* Q (i,f),* t/3 (t,f )} =8(z- z')X a (r,6)xW{r',6'), (58) 
{X a (t,x),X^(t,x')} =5(z-z')Y a (r,e)Y^(r',e'), (59) 



(57) 
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with a and /3 the spinorial indices, and all the other anticommutators vanishing. With equation ( p5| ) and these 
anticommutation rules, it follows immediately, for creation and annihilation operators, that 

{b{-k),tf(-k')} = {b{-k),tf(-k')\ = 27r2fc2fc , ( 5(fc - A;'), 

(60) 

[d[k),ctf{k')\ = {d(k),<tf(k')} = 2Tr2k2k'5(k-k'), 

with all other anticommutators vanishing. From the expressions ( |50"| ) and ( pi] ) and with the anticommutation rules 
(|60|), it is possible to derive the anticommutator between two quantum field operators at any time. For instance, the 
anticommutator between 'J and W reads 

{*„(*), * f V)} = / (2 ^22fc2fc^ { bt (" fc )^ (t+Z) +b(-k)e-^\b(~k')e-^^' + ^ + 

+ 6 t (-fc')e <k ' (t ' + * /) } Jra^fljX^^.fl 7 ). (61) 

Thanks to the delta function coming from the anticommutators between the b and w , this expression reduces to 

{* a (x), *t^(x')} = id t A(t -t' + z- z')X a (r, 6)XW(r', 9'), (62) 

with A(x) the well known Pauli- Jordan function which vanishes for spacelike separation, so the spinor fields indeed 
respect micro-causality along the string. 



4- Fock states 



The Fock space can be built by application of the creation operators on the vacuum state |0) which by definition 
has to satisfy 

6(-*)|0) = 6(-fc)|0) = d(k)\<D) = d(k)\%) = 0, (63) 

and is normalized to unity, i.e., (0|0) = 1. Each Fock state represents one possible combination of the fields exitation 
levels. Let \V) be a Fock state representing particles labeled by i and N$ anti-particles labeled by j, of kind ^, 
with respective momenta fcj and lj, and, N x particles labeled by p and N x anti-particles labeled by q, of kind X, with 
respective momenta r p and s q . By construction the state is 

\V) = b\-h) . ..b\-k t ) . . . bH-k N Jb_H-h) ■ • .6 + Hi) • ..b\~l w Jd^ n ) . . . 

. ..d\r p ) . ..d\r Nx )S( Sl ) . ..d}(s q ) . . . d f (3^)10). (64) 

Normalizing such a state is done thanks to the anticommutators (|60|). For instance, for a one particle 'J state with k 
momentum, using Eq. (|63|), the orthonormalization of the corresponding states reads 

(k'\k) = 2w2k2k'S(k- k'). (65) 

Obviously similar relations apply to all the other particle and anti-particle states. Keeping in mind that the observable 
values of quantum operators are their eigenvalues in a given quantum state, let us compute the average of the 
occupation number operator involved in many quantum operators, as will be shown. For ^ particles it is 

and analogous relations for the other particle and anti-particle states. By definition of the Fourier transform (^2h, the 
infinite factor 6(0) is simply an artifact related to the length of the string L by 

L = 27n5(0), (67) 

in the limit where this string length L — > oo. Note that using periodic boundary conditions on L allows to consider 
large loops with negligible radius of curvature. 
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B. Fermionic energy momentum tensor 



The simplest way to derive an energy momentum tensor already symmetrized is basically from the variation of the 
action with respect to the metric. Moreover, the Hamiltonian density of the fermion "J, 

H^=d t Il^ + d t m il -£^, (68) 

with IT the conjugate field II = 1^7°, is also equal to the T" component of the stress tensor. In our case the metric 
is cylindrical and we assume a flat Minkowski space-time background, thus in the fermionic sector the stress tensor 
reads 



Tr =2-^-g^d > and = 2—^ - g^C x . (69) 

v og^v * og^ 

Once again, let us focus on 'J'. Plugging Eq. (||) into Eq. (|6^) gives 

I*" = 1*7^9^* - - (<9 ( ^) 7 l/) * - B^fJ. (70) 



1. Symmetrized Hamiltonian 



From Noether theorem, the Hamiltonian P* is also given by the conserved charge associated with the time component 
of the energy momentum tensor 



T" = i^dt^ - i (<9 t ¥) 7 *. 



(71) 



Thanks to the expression of the quantum fields in equations (|50| ) and ( |5l|) , and using the properties of the zero modes 
from Eq. (45), the quantum operator associated to Tl* reads 



Titt 
1/1 



dfcdfc' 
{2n) 2 2k 

+ (tf (-k)b(-k') - &(-fc') &t (- fc )) e" i(fe ' 



(-b{-kp{-k') + bH-k')b(-k)^j e w-m+z) 

-i(k' -k)(t+z) 

+ (b{-k)b{-k')-b{-k')b{-k))c- l{k ' +k){t+z) 



+ (-b\-kp(-k') + b^(-k')b^-kfje^ k ' +k ^ t+z ^ \X\ 



The Hamiltonian is given by spatial integration of the Hamiltonian density, or similarly from Eq. ( |72| ) , 

dk 



P, 



2n2k 



-b(-k)tf(-k)+tf{-k)b(-k) 



(72) 



(73) 



Note, once again, that all the terms in the form 6^6+ or bb vanish as a consequence of the chiral nature of the spinor 
fields which only allows k > 0. The average value of this Hamiltonian in the vacuum is not at all positive, but a 
simple renormalization shift is sufficient to produce a reasonable Hamiltonian provided one uses fermionic creation 
and annihilation operators with the corresponding definition for the normal ordered product (antisymmetric form). 
The normal ordered Hamiltonian is therefore well behaved and reads 



dfc 
27r2fc 



6 + (— fc) + b}(-k)b(-k) 



(74) 



However, note that such a normal ordering prescription overlooks the differences between the vacuum energy of empty 
space, and that in the presence of the string for the massless fermions. Formally, from Eq. (|72|), the normal ordered 
Hamiltonian is also obtained by adding the operator J da? V,/, to the infinite Hamiltonian in Eq. (|73]), with 



Vj, =: T, h 



Tit 



dfcdfc' 
(27r) 2 2fc 



{b{~k),b\~k')}S k -*)(*+-' 



+ {&(-fc'),6 t (-fc)}e- l ( fe '- fc )(*+^j |X| 5 



(75) 
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Owing to the anticommutation rules in Eq. (pq), this expression reduces to 

dfc 



Vv, = / — k\X\ z . (76) 



7T 



This infinite renormalizing term of the vacuum associated with the zero modes on the string comes from the contri- 
bution of the infinite renormalization of the usual empty space together with a finite term representing the difference 
between the two kinds of vacua. The previous expression (|76| ) emphasizes the structure of the divergence, and it 
can be conjectured that the finite part is simply obtained by a cut-off A in momentum values. The finite vacuum 



contribution to the stress tensor can therefore be represented from Eq. (76), up to the sign, by the energy density 



^\X\' = ^\X\\ with L v = ^. (77) 

The precise determination of the value of L v is outside the scope of this simple model. It is well known however, that 
the vacuum effects generally involve energies smaller than the first quantum energy level and consequently it seems 
reasonable to consider that L v > L. For a large loop, L v can be roughly estimated using the discretization of the 
momentum values. With k n — 2irn/L, therefore reads 

A 00 
n=0 

Assuming that the vacuum associated with the fermionic zero modes on the string matches the Minkowski one 



associated with massless fermionic modes, in the infinite string limit |31 3^ , L v can be obtained by substracting the 
two respective values of Vy,, once the transverse coordinates have been integrated over. The infinite sum over n can 
be regularized by a cut-off factor e _efc ™ , letting e equal to zero at the end of the calculation |3j| . The regularized 
expression of finally reads 

- oc 

Vv- = ^l* 2 |]>> e " £ *"> (79) 
and expanded asymptotically around e — 0, it yields, once the transverse coordinates have been integrated, 



rdrd*V,~^-^. (80) 



As a result, the infinite renormalizing term relevant with the usual vacuum associated with two dimensional chiral 
waves is just 1 / ire 2 whereas the relevant vacuum associated with the zero modes along the string is exactly renormalized 
by J rdrd9V^ given in Eq. (|30|), and therefore involves the finite term n/3L 2 with a minus sign. As a result, the 
string zero mode vacuum appears as an exited state in the Minkowski vacuum associated with two dimensional chiral 
modes, with positive energy density ir/3L 2 . In this case, the short distance cut-off therefore reads 

L 2 = 6L 2 . (81) 

It is therefore necessary to add the 2irL / L 2 term to the canonical normal ordered prescription, in Eq. ([74j) , in order to 
obtained an Hamiltonian with zero energy for the Minkowski vacuum associated with two dimensional chiral modes. 

It is important to note that, in order to be consistent, the previous calculations can only involve the vacuum 
associated with the corresponding quantized modes, i.e. in this case the zero modes. Therefore this does not take 
care of the massive modes. In fact, even for zero occupancy of the massive bound states, the physical vacuum along 
the string must involve a similar dependence in the vacuum associated with the massive modes. The influence of the 



two dimensional massive vacuum on the equation of state will be more discussed in section VB 



2. Stress tensor 



All the other terms of the energy momentum tensor can be derived from equation (]69|). From the relationships 
verified by the zero modes currents in Eq. (Eq), all the transverse kinetic terms vanish. Moreover, the only non- 



vanishing components of the axial and vectorial currents (27) are A and Finally, the energy momentum tensor 
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reads 





rpte 
\ T tz 



T$ 


Tf 







T, 



1> 



(82) 



rpZZ I 
1 li> I 



Because the zero modes are eigenvectors of 7°7 3 , the operators 7°9° and 7 3 <9 3 are formally identical for each spinor 
held, and therefore the diagonal terms of the stress tensor are identical 



T'ZZ 



r'tt 



(83) 



From equation (|69|), we find the transverse terms to be T™ = — B 
components yield the transverse terms T 



rpty 



9 At 



and Tf = —B ( 



In a Cartesian basis, these 



rpZX 



and , which vanish once the transverse degrees of freedom 



have been integrated over. The only non-vanishing non-diagonal part of the stress tensor comes from the lightlikc 
nature of each fermion current and reads 



(84) 



while the counterpart of the field X gets a minus sign because of its propagation in the 

Tl z = 2iXj d f X. 



-direction. 



(85) 



In the above expressions, the backreaction is neglected, but the fermionic current, j** = + j£, generates, from 
Eq. (|2l|), new gauge field components, B t and B z , which have to be small, compared to the orthoradial component, 
B@, in order to avoid significant change in the vortex background. However, up to first order, they provide corrections 
to the energy momentum tensor whose effects on energy per unit length and tension will be detailed, in the classical 
limit, in section [v|. The backreaction correction to the two-dimensional stress tensor then reads 



Iff = 



1 -2B t f 



1 



2 
Btj 



{d r B t f 



; (d r B z ) 2 B z f - B t f - (d r B t )(d r B z ) 



{d r B t )(d r B z ) 2B z f 



~ 2 idrB z f 



(86) 



C. Axial and vectorial currents 



From the expression of the spinor fields, the current operators are immediately found in the Fock space. Moreover, it 
is interesting to compute the electromagnetic-like fermionic current, its scalar analogue being involved in the equation 
of state for a cosmic string with bosonic current carriers Jl4] |. From an additional global U(l) invariance of the 
Lagrangian, the electromagnetic-like current takes similar form as the vectorial one coupled to the string gauge field. 
It physically represents the neutral limit of the full electromagnetic coupling. 



1. Vectorial currents 



Let J* 1 be the electromagnetic current in the neutral limit. From Noether theorem with global U{1) invariance, 
this is 

= + J£ = -* 7 ^* + Xj^X, (87) 

The fermions 4" and X carry opposite electromagnetic- like charges in order to cancel anomalies pi. Using equation 
(|50|), their components read 

■4 : = - : J$ ■■ = ■■ J4, ■■ \X\ 2 , 



J x : — : J x : — : J x : \Y\ 



with the quantum operators defined as 
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-b\-k')b(-k)e^ k '^ t+z ^ +bH-k)b{-k')e~ l{k '- k){t+z) + 



dkdk' 
(27r) 2 2fc2fc' 

+ b{-k)b{-k')e-^ k+k '^ t+ ^ +6 t (-fc)fe t (-fc')e i(fe+fc ' )(t+z) j \X\ 2 , 

= / (2^2*2* [-^(feO^e^'-^^^+^CWOe^^'^' 
+ rf(fc)d(fc')e- l(fc+fc ' )( *- 2) +rf t (-fc)d t (fc')e l(fc+fe ' )( *- 2) l |F| 2 . 



(89) 



(90) 



The conserved charges carried by these currents are basically derived from spatial integration of the corresponding 
current densities, and are 



dk 



2ir(2k) 2 
dk 



27r(2fc) 2 



b\-k)b{-k) + ti}(-k)b(-k) 
-d f (k)d(k) + d\k)d(k) 



(91) 
(92) 



As expected at the quantum level, the anti-particles carry charges that are opposite to that of the particles for both 
fields. This is again because the opposite chirality of the fields makes the bb and ww terms vanishing. The vectorial 
gauge currents are easily obtained from the neutral limit ones replacing the electromagnetic charge by the U(l) gauge 
one, namely 



'Xv 



q C ^+ C ^ :J^:\X\\ 
q C *»+ C ^ : J y : \Y\ 2 . 



(93) 



Axial currents 



In the same way, the axial currents are derived from their classical expressions as function of the quantum fields, 
from Eq. j45|), 



t 



J 



XA 



Hi 



J XA 



q 



-Xh ^Xl 



■■Jx--(\ti\ 2 -M 2 )- 



(94) 



Thanks to the normalizable zero modes in the transverse plane of the string, it is possible to construct a Fock space 
along the string. The chirality of each spinor field being well defined, anti-particle states appear at quantum level as 
another mode propagating at the speed of light in the same direction than the particle mode, but carrying opposite 
gauge and electromagnetic-like charges. The observable values of the quantum operators previously defined are given 
by their average value in the corresponding Fock state. In particular, the energy per unit length, the tension, and the 
current per unit length, can now be derived from the previous expressions. 



V. EQUATION OF STATE 

In the case of a scalar condensate in a cosmic string, it was shown by Peter |l4| that the classical formalism of 
Carter (l9) with one single state parameter could apply and an equation of state for the bosonic cosmic string could 
be derived in the form 



U-T = s/\w\C, (95) 

where U and T are respectively the energy per unit length and the tension of the string, C is the current density along 
the string and w a state parameter which appear as the conjugate parameter of C by a Legendre transformation. 
An analogous relation can be sought for our string with fermionic current-carriers from the classical energy per unit 
length, tension, and current density values. 



Consider the fermionic cosmic string in the quantum state (64). The energy per unit length and tension in this 



state are basically given by the eigenvalues associated with timelike and spacelike eigenvectors, respectively, of the 
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average value in \P) of the energy momentum tensor, once the transverse coordinates have been integrated over. The 
stress tensor is obviously the total energy momentum tensor 

T"" = T% v + Tj?+ : T$ v : + : Tg v :, (96) 

where Tg v and T£ v are the gauge and Higgs contributions which describe the Goto-Nambu string and which integrated 
over the transverse plane provides only two opposite non- vanishing terms 

T dr d6 (Tf +T%) = - J r dr d6 (T 22 + T h 22 ) = M 2 , (97) 

thus defining the unit of mass M. 

A. Average values in the Fock state \V) 

1. Two-dimensional energy momentum tensor 

— a/3 

Now, let us define the energy momentum tensor operator in two dimensions, T say, once the transverse coordinates 
have been integrated over, and where we have suppressed the corresponding vanishing terms. Therefore, with a and 
(3 equal to t or z, and neglecting for the moment the backreaction, it reads 



T 



: „ , , M 2 + Jr drd9 (: T« : + : T» :) Jr dr d9 (: T^ z : + : 



Jr drdO ( : T^ z : + : T l x z :) —M 2 + frdr d6 (: Tff : + : T** 



The average value in the Fock state \P) of , is immediately obtained from equations (66) and ([72| 



-a/3 



^a/3 N (V\T H \V) _ ( M 2 + E Xv + E^ E Xv - E in 



[ )V (V\V) \ E xv -E^ v -M 2 + E xv +E^ I' 



with the notations 



E xv = I r dr d6{: T« =>7> = 7 ( J> + £ l i ] > ( 10 °) 



'jv v 




, P =i 



The summations are just over the momentum values taken in each particle and anti-particle exitation states, so that 
E Xv and E^ v depend on the quantum state \V). The 2ttS(0) factor has been replaced by the physical length L in 
order to deal only with finite quantities. Moreover, from the integral expression of the normal ordering prescription 
in Eq. (|77|), the quantum zero mode vacuum effects appear simply as a shift of the previous expressions, and the 
corrected values of the parameters E^ v and E Xv therefore read 

2ir _ 2n 

Eip-P v = EiPv + "72 ' ano - ^xp, = E Xv + (102) 

l-i., L,, 



Note that, from equations (10C), (101) and (102), if L v > L, the vacuum contribution can be neglected for non-zero 
exitation states. 

2. Current densities 

In the same way, the average value of current operators in the Fock state \V) are basically derived from the average 
of the operators : : and : J x : 

(: Ji, :) V = -N i , + N. 4 „ (: J x :) v = -N x + N x . (103) 
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Therefore, the electromagnetic-like current per unit length in the neutral limit becomes, after transverse integration, 

(JV - 1 [{-N^ + N^) + (N x - N x )] , (104) 

(T) v = i [(^ - 7v^) + (tv x _ ]v x )] . (105) 



Averaging the vectorial and axial gauge currents in equations (|93j) and (|9J) allows a derivation of the total fermionic 
gauge current density 

<j*)v = \ [Mr) (-N^ + N^) + f x (r) (-N x + N x )] , (106) 

(f)v = \ [Mr) (N^ - N*) + f x (r) (-N x + N x )] , (107) 

with the radial functions 

fx(r) = 9c XR \^\ 2 + qc XL \U\ 2 , (108) 

Mr)=qc^ 2 \ 2 +qc^ 3 \ 2 . (109) 

Moreover, note that these currents can be lightlike, spacelike or timelike according to the number of each particle 
species trapped in the string. For instance, the square magnitude of the electromagnetic-like line density current reads 

C% = (f)% - ( T% = (jj (N^N X + N^N X - N^N X - N X N^) . (110) 

As expected, if there is only one kind of fermion, 4" or X, which respectively means N x = N x = or N$ = Nj, = 0, 
the current is lightlike. However spacelike currents are also allowed from the existence of anti-particles as they result 
from simultaneous exitations between particles of one kind and anti-particles of the other kind (as for instance ^ 
and N x ^ 0, or N x ^ and ^ 0). Finally, timelike currents are obtained from simultaneous exitation between 
particles or anti-particles of both kind (N x ^ and N$ ^ 0, or N x ^ and ^ 0). 

B. Energy per unit length and tension 

In the case of a string having a finite length L, periodic boundary conditions on spinor fields impose the discretization 
of the momentum exitation values 

h - 2lT i - 2lT - - 2n - 2n - run 



wher e ru )i , •n n j )j , n Xp and n Xq are positive integers given by the particular choice of a Fock state. From Eq. ( Jl 00| ) and 



Eq. (101), the parameters E Xv and E^ v therefore read 





and E^ v =7? V n * + V n * ■ (H 2 ) 



In the preferred frame where the two-dimensional energy momentum tensor is diagonal, the energy per unit length 
and the tension appear as the eigenvalues associated with the timelike and spacelike eigenvectors, respectively. By 
means of equation (199) , they read 



U-p = M 2 + 2^/E~E~~, (113) 

T P = M 2 - 2y/E Xv E^. (114) 
Note, first that the line energy density and the tension always verify (22) 

U + T = 2M 2 . (115) 



Moreover the zero mode vacuum effects just modify the parameters E Xv and E^ v as in Eq. (|102| ), and therefore 
do not modify this relationship. On the other hand, massive modes, because they are not eigenstates of the 7°7 3 
operator, yield vacuum effects which certainly do not modify the time and space part of the stress tensor in the same 
way, as was the case for the zero mo des [ see Eq. (p3|)]- As a result, it is reasonable to assume that the massive mode 
vacuum effects modify the equation flll£| ) by just shifting the right hand side by a finite amount, of the order 1/L 2 . 
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1. Classical limit for excited strings 



In order to derive classical values for the energy per unit length and tension, we do not want to specify in what 
exitation quantum states the system is. If the string is in thermal equilibrium with the external medium, it is 
necessary to perform quantum statistics. The number of accessible states in the stri ng is precis ely t he total number of 



combinations between the integer n^, n^-, n Xp , and n x , which satisfies equations (IOC) and (101) for fixed values of 
the stress tensor, or, similarly, at given E x and E^. A possible representation of such equilibrium is naturally through 
the microcanonical entropy 

S = k b lnQ, (116) 

with D the number of accessible states and fcf, the Boltzmann constant. Let Q(p, N) be the well known partition 
function Q which gives the number of partitions of the integer p into exactly N distinct non-zero integers. With the 
following integers 



K x = ^ E x> and K * = ^ E 4- (ll7! 



the number of accessible states f2 reads 

l? x (~ x + l) 

JV V n 



n= E Q(rbc> N x) E Qfcx> N x) 



_ N X (N X + 1) 
* 1 



Ktjj ^ K-^—riii 



E 0(n^,JV«) E Qin^N^,). (118) 



The energy per unit length and tension of the string will therefore be the values of U and T which maximize the 
entropy at given , , N x and N x . This formalism might be useful whenever one wants to investigate the dynamics 
of the string when the massless current forms, i.e., near the phase transition at high temperatures. In what follows, 
we shall assume that whatever the mechanism through which the fermions got trapped in the string, they had enough 
time to reach an equilibrium state with vanishing temperature. This can be due for instance by a small effective 
coupling with the electromagnetic field opening the possibility of radiative decay p2| . If no such effect is present, 
then one might argue that the string is frozen in an exited state, the temperature of which possibly playing the role 
of a state parameter jl9| for a macroscopic description |36j . 

Note that for a given distribution such as those we will be considering later, the occupation numbers at zero 
temperature must be such that, owing to Pauli exclusion principle, the interaction terms implying for instance a \D r \f r 
decay into a pair XX through Higgs of exchange are forbidden (vanishing cross-section due to lack of phase space). 
In practice, this means that the following analysis is meaningful at least up to one loop order. 



2. String at zero temperature 



For weak coupling between fermions trapped in the string and external fields, as is to be expected far below the 
energy scale where the string was formed, the set of particles is assumed to fall in the ground state and because of 
anticommutation rules (60) it obeys Fermi-Dirac statistic at zero temperature. Consequently, the parameters E Xv 
and E^ v reads 



2ir _ 2tt 

E i> = ~L [P^( L Pi> + !) + P^,i L P^ + !)] > and E x = Y \-Px( L Px + !) + P x ( L Px + 1 )] ' 



(119) 



where the new parameters p = N/L are the line number densities of the corresponding particles and anti-particles 
trapped in the string. Strictly speaking, the se ar e the four ind ependent state parameters which fully determine the 
energy per unit length and the tension in Eq. ( |l!3| ) and Eq. (114), and so cosmic strings with fermionic current-carriers 
do not verify the same equation of state as the bosonic current-carrier case. This is all the more so manifest with 
another more intuitive set of state parameters, 1Z and O, defined for each fermion by 
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These are simply polar coordinates in the two-dimensional space defined by the particle and anti-particle densities 
of each spinor field. The parameter 1Z physically represents the fermion density trapped in the string regardless of 
the particle or anti-particle nature of the current carriers. It is therefore the parameter that we would expect to 
be relevant in a purely classical approach. On the other hand, O quantifies the asymmetry between particles and 
anti-particles since 



p = 7\Lcos0 — — , and p = 1Z sin 9 — 
The energy per unit length, tension and line density current now read 



U = M 



T = M 2 




C 2 = 87^7^ sin (q x - ^ sin (o 



(121) 

(122) 

(123) 
(124) 



There are always four independent state parameters but only two, 7Z X and TZ^, are relevant for line density energy 
and tension. Compared to the scalar case where only one kind of charge carrier propagates along the string, it is not 
surprising that we found two degrees of freedom with two kinds of charge carriers. On the other hand, the nature of 
the line density current is not relevant because it only appears through 0^, and O x , which not modify U and T, at 
least at the zeroth order. The energy per unit length and tension relative to M 2 are represented in Fig. [j] and Fig. || 
as function of 7Z X /M and TZ^/M, in the infinite string limit. 




0.5 



FIG. 1. The energy per unit length, in unit of M and at zeroth order, as function of X and fermion densities plotted in 
unit of M, in the infinite string limit. 
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0.5 



FIG. 2. The tension, in unit of M 2 and at zeroth order, as function of X and fermion densities plotted in unit of M, in 
the infinite string limit. 

As expected from their analytical expressions in the infinite string limit, the energy per unit length is always positive 
and grows with both fermion densities, 1Z X and TZ^, whereas the tension always decreases and takes negatives values 
for large fermion densities. Obviously, in the case 1Z X = 7Z^ = there is no current along the string and we recover 
the Goto-Nambu case, U = T = M 2 . The chiral case, where the fermionic current is lightlike, is o btained for 1Z X = 0, 



or 1Z^ = 0, and also verifies U = T = M 2 as in the chiral scalar current case pq] . From Eq. (123), and in the infinite 
string limit, the densities for which the tension vanishes verify 

M 2 

n xK<p = (125) 

This curve separates the plane (7Z X , TZ^) in two regions where T is positive near the origin, and negative on the other 
side (see Fig. ||). In the macroscopic formalism of Carter fLO] , the transverse perturbations propagation speed is given 
by dp = T/U, and therefore the domains where T < correspond to strings which are always locally unstable with 
respect to transverse perturbations. 
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FIG. 3. Sign of the tension in the (1Z X /M, IZ^, /M) plane, in the infinite string limit. According to the macroscopic formalism, 
the string is unstable with respect to transverse perturbations for T < 0. 

The tension of the string becomes negative only for carrier densities close to the mass of the Goto-Nambu string 
M. For such currents, it is necessary to derive the backreaction in order to see how relevant it is for the energy per 
unit length and tension. Moreover, in a renormalizable model, the vacuum mass acquired by the fermions from their 
coupling to the Higgs field is less than the Goto-Nambu string mass, and thus, another quantum effects may take 
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place before the negative tension is reached, like tunneling into massive states. Besides, note that M, the string unit 
of mass, arising from non-perturbative effects, may well be much larger than the Higgs mass, and so, it is expected 
that "R X ,1Z^ <C M. Thus, the no-spring conjecture Q proposed in the case of bosonic carrier presumably apply to 
the fermionic carrier cas e as well. Moreover, the zero mode vacuum effects on energy per unit length and tension 
appear clearly from Eq. (102) as additional string length. The corr ected values of the equation of state are therefore 
obtained by replacing the physical length of the string, L, in Eq. (122) and Eq. (123), by an equivalent length, L e 
say, which verifies 



1 

2Z! 



1 

2L 2 



1 

a' 



(126) 



In the particular case where L 2 = 6L 2 , it reads L 2 = (3/2)L 2 . On the other hand, the massive vacuum effects 
certainly shift in a different way U and T by a finite amount as previously discussed, but will not be considered in 
the following. In the next section, the backreaction is derived in the classical limit in order to find corrected values 
of energy per unit length and tension. Moreover we shall take care of the finite length of the string L, keeping in 
mind that its value, and consequently the value of L e , have to be larger than X/M since all physical values have been 
derived in the classical vortex background, i.e., the quantum effects of the Higgs field have been neglected. 



C. Backreaction 



The existence of fermionic currents carryin g ga uge charge along the string gives rise to new gauge field components, 
B t and B z , from the equations of motion (plf). These, being coupled with the corresponding currents, provide 
additional terms in the energy momentum tensor (pq). As a first step, the new gauge field components are computed 
numerically from the zero modes solutions of Eq. Q32| ) . The corrected equation of state is then analytically derived, 
the numerical dependencies having been isolated in model dependent coefficients. 



1. Backreacted gauge fields 

In order to compute the B t and B z fields at first order, we only need the zeroth order values of the zero modes and 
the vortex background. Let us introduce the dimensionless scaled fields and variables 

<p = r}H, Qe = Q and r = — , (127) 

m h 

with rrih = rjy/X the classical mass of the Higgs field. From the equation of motion (|l6|), the orthoradial gauge field 
Q and H are solution of 

\H{H 2 -1), (128) 



d 2 H 


1 dH 


HQ 


dg 2 4 


g dg 


B 2 


d 2 Q 


1 dQ 


= m| j 


dg 2 


g dg 


m h 



>, (129) 

where mj, = qc§r\ is the classical mass of the gauge boson. The numerical solutions of these equations have been 
computed earlier by many people using relaxation methods |29|. They are presented in Fig. |^ for a specific 

(assumed generic) set of parameters. 
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FIG. 4. The solutions of the field equations for the vortex background. The Higgs field, H, takes its vacuum expectation 
value at infinity and the gauge bosons condensate in the vortex. 



In the same way, deriving Eq. 



with respect to g yields the right component of the zero mode Y as a solution of 



Cl 1 H + e 



C0 



Q-nH' Q 



Q 



H 



C Xl c XR ( Q n 



rn i 



9L 

Q 



- -4h 2 Ci = o, 

m h 



while the left one satisfies 



i mf TJt -t> C * R 



-Ci, 



(130) 



(131) 



where a prime indicates a derivation with respect to the dimensionless radial variable g. The field ^> verifies similar 
equations with the transformation, £ — * £ and c x — > c^. The numerical integration has been performed with a 
relaxation method |^9| and verified on the original system ([52]) with a shooting method. As a result, the normalized 
probability densities of the zero modes, \X\ 2 and \Y\ 2 , are plotted in Fig. [|. The dimensionless radial functions f x 
and defined from Eq. (108) and Eq. (109) by 



fx 



2tt 



■f: 



2 JXi 



and — — J ftp j 
m h 



(132) 



are plotted in Fig. ^. As expected, the fields are confined in the string core, and so will the corresponding fermionic 
currents. 

Let us define the more relevant components of the backreacted gauge field, AB = B z — B t and T,B = B z + B t , with 
the corresponding dimensionless scaled fields AQ and EQ defined by 



ml AQ 



AB = (p - p x ) 



and EB 

The equations of motion ( pl| ) in the classical limit now reads 



W — I ■ 



f: 



AQ" + -AQ' -^H 2 AQ = ^, and EQ" + -EQ' - ^|iJ 2 EQ = 



'"7, 



9C0 



'"7 



(133) 



(134) 



As for fermions, these new gauge fields get their masses from coupling with the Higgs field, and therefore have non- 
zero mass outside the string core. Moreover, they are generated by fermionic massless currents confined in the core , 
therefore they also condense in and do not lead to new long-range effects. The solutions of these equations (134) 
have been obtained using, once again, a relaxation method |29| and are represented in Fig. Note that owing to 
the scaled field AQ and EQ, we have separated the numerical dependence in the gauge field and currents from the 
fermion densities content [see Eq. (|133|)]. 
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FIG. 5. The normalized probability densities of the zero modes, \X\ 2 and \Y\ 2 . The rapid decay far from the string core 
reflects the bound state nature of the condensates. 
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FIG. 6. The dimensionless radial current functions f x and fy. They can be viewed as the effective transverse density charge 
carried by the fermion currents. Their sign results in the initial choice of each conserved fermion gauge charge. 
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FIG. 7. The dimensionless backreacted gauge fields AQ and SQ, generated by the fermion currents. They do not lead to 
new long-range effects since they acquire non-zero mass outside the string core due to their coupling with the Higgs fields. 
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Up to now, we have computed the fermionic gauge currents along the string as well as the component B t and B z , 
so that the backreaction correction to the energy momentum tensor is computable from Eq. . As in the previous 
section, the energy per unit length and tension can be derived in the preferred frame where the stress tensor is 

diagonal, but now we have to find the eigenvalues of the full two-dimensional energy momentum tensor T b r + (T ), 
with 

T£j. = f rdrdOT^ (135) 



2. Energy per unit length and tension with backreaction 



Using the dimensionless field, AQ and T,Q, with the expressions of the currents given in equations (106) and ( 107 ) 
one gets, after some algebra, the full expression of the stress tensor with corresponding eigenvalues 



u = m 2 - i(Q X) &^n x % + 4tt J (n x _ j _ _1_ ] , , m 



f=M--HU,.U,.YR x R„ Jm/(^-2I|) (^-2lf)' (137) 



with the scaled state parameters 



TZ X =U x Jl- ^ (I Ax Ia2) sin 2 (9 X - -) , (138) 



1Z 4 , = K^Jl- ^2 (J Ev > - 7 S2 ) sin 2 [e t - -J , (139) 

and the function I{Q X , @ip) defined by 

/( } = ml (7 Sx+ / A ,)sin(e x - f )sin(e,- f) 

7T77 1 ' i .... 



1 - ^ ^ - IA2) sin 2 (0 X -1)^/1-^ (J EV , - / S2 ) sin 2 (0* - f 

The numerical integrations previously carried out appear through pure numbers which depend only on the model 
parameters. The coupling B^j^ leads to the following quantities 

J Sx = -2 / gdgXQig)^, I M , - -2 / gdgAQ(g)^, (141) 



9C0 J qc^ 

I^ = -2 I gdg^Qig)^^, I A% = -2 / gdg AQ(g)^, (142) 



qccj, J qc^ 

while the kinetic contribution of the new gauge fields appears through 



7 A2 = J gdg (d e AQ{g)) 2 , 7 S2 = / gdg (d e £Q(g)) 2 . (143) 



By means of the equations of motion (134) and the constant sign of AQ and T.Q, I& x and Ij^, are found to be always 
positive. Intuitively, as in electromagnetism, the gauge field generated from charge currents tends to resist to the 
currents which give birth to it. In our case, the backreaction actually damps the weight of the charge carriers in the 
energy per unit length and tension. In fact, the relevant state parameters are now 1Z instead of 1Z with 1Z < TZ s ince 
Ia x and are positive. Moreover, numerical calculations show that the kinetic contribution numbers ( |l43|) ar e 
always one order of magnitude smaller than those resulting in the coupling between gauge fields and currents ( |141| ) . 
as expected for reasonable backreacted gauge field since they only involve the square gradient of these fields [see 



Eq. (143)]. However, there is an additional term involving new dependence in the asymmetry between particles and 
anti-particles through the I(O x , Qf) function. In order to understand this point physically, let us derive the magnitude 
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of the gauge current carried by the fermions. From Eq. fllOq ) and Eq .(107), once the transverse coordinates have 
been integrated over, the dimensionless magnitude reads 



^2 

i 



2(2g C0 ) 2 



F X F^ sin (8 X - f ) sin (6^ - § )U X K V 



v±* - ^ sin2 ( e * -f)\A 

with the dimensionless constants are 
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These numbers can be viewed as the effective charge carried by the fermionic gauge currents since 

T 2 



^2 
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(2qc^fF x F^ 
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T 



The function I(O x ,0,/,) therefore verifies 



8nF v F 



I Alb ^ 2 



(144) 



(145) 



(146) 



(147) 



and, as before, according to Eq. ( |134| ), {Iy, x + lAip)/{F x F^) is always positive, so the sign of / directly reflects the 
spacelike or timelike nature of the current. Thus, in addition to the backreaction damping effect, there is a correction 
to the energy per unit length and tension directly proportional to the magnitude of the fermionic current. Note that 
this effect appears as a correction due to backreaction and not, as it is the case for cosmic string with bosonic current 
carriers, at the zeroth order EM. 



3. Equation of state with backreaction 

Unfortunately, the corrected expressions of the line density energy and tension involve four independent state 
parameters, and consequently are not easily representable. However, they can be studied as functions of the damped 
fermion densities 1Z, modified only by the function I{Q X , O^,) which quantifies the efficiency of the fermionic currents 
in generating backreactcd gauge fields. In this way, the comparison with the zeroth order case is all the more so easy. 

The study of the surfaces defined by U and T in the plane (TZ Xl 'R-^) is less canonical than at zeroth order. Three 
critical values of the function / are found to modify the behaviors of the tension and energy per unit length, namely, 
— 4-7T, 0, and 4-7T. However, only small values of / are reasonable in this model as it is discussed in the next section. 
This analysis is consequently constraints to values of |/| < 4ir. 

a. Energy per unit length. The line density energy follows different behaviors according to the value of J. 

The first and simplest case / < 0, obtained for spacelike fermionic gauge currents, is very similar to the zeroth 
order case, and the energy per unit length just grows a bit faster with the damped fermion densities 1Z X and TZ,p, as 
on Fig. ||. 

For timelikc currents, / > 0, we find that the backreaction damps the growth of the density line energy with the 
fermion densities. As a result the line density energy seems to decrease in some regions, and the stationary curves of 



U with respect to TZ X are given, from Eq. (136), by 

dU 




0^7^ = — Tl Xi — 2 — - — s- , (148) 



dK x L e *y + (i 6 ^2 _ p )U 2 



and thanks to the symmetry between 1Z X and IZ^p, similar equations are obtained for dU/dlZ^, = 0. Finally, the 
variation domains of the line density energy are represented in Fig. [)] for < / < 47r. The first discrete values of 
the fermion densities (the length of the string is finite) have been represented by dots in the (7?. x ,7?.^)-plane, and 
as can be seen in Fig. |9|, for reasonable values of /, there is no available quantum state inside the tiny decreasing 
regions. Consequently, the density line energy always grows with the fermions densities and remains positive. Since 



the stationary curves of U are asymptotically proportional to 1/L e < l/L [see Eq. (148)], they coincide with the axis 
in the infinite string limit. The surface describing U(1Z X ,1Z^) has also been plotted in Fig. ^| in unit normalized to 
M 2 , and for minimal acceptable value of L e = 10/M just in order to show the influence of the finite length. 



24 



0.5 

FIG. 8. The energy per unit length in unit of M 2 as function of 1Z X /M and IZ^, /M, for spacelike currents with /(© X5 ©v) < 0- 
The influence of the finite length of the string just appears near the axes. 
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FIG. 9. Variation domains of the energy per unit length, plotted in unit of M 2 , in the (iZ x /M,1Z^/M) plane for timelike 
currents with < I < 4tt. The vertical hatched regions define domains where the line density energy could decrease with 1Z X 
whereas the horizontal ones are regions where the line density energy could decrease with TZ^,. On one hand, these domains 
are asymptotically limited by the curves 1Z = (2ir/L e ) -y/2/(167r 2 — I 2 ) and therefore coincide with the axis for large values of 
L e - On the other hand, for reasonable values of / 4n, there is no accessible quantum state inside, the first one being shown 
as a dot. As a result the energy per unit length always grows with the parameters 1Z and is always positive. 



b. Tension. The study of the tension with respect to the fermion densities is performed in the same way. As 
before the stationary curves of T with respect to 1 Z V or IZ^p are found from Eq. (137), and follow the same equation 
as those of the energy per unit length in Eq. (148), although the variation domains are not the same and have been 
plotted in Fig. |ll| for different values of the function /. 

For timelike fcrmionic gauge current, I > 0, the tension decreases faster than in the zeroth order case, with the 
damped fermionic densities 1Z X and H^, as on Fig. 10, and reaches negative values at large densities (see Fig. [l2"| ). 
The backreaction just increases the slope of the surface, and thus, the negative values are reached more rapidly. As 
for the energy per unit length, the equivalent length was chosen equal to L e — 10/M in the following figures. 
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FIG. 10. The tension for timelike currents with / > 0, plotted in unit of M 2 in the (lZ x /M,lZ,p/M) plane. Note the shift 
near the axes due to the finite length of the string. 

For spacelike fermionic gauge currents, — 47r < I < 0, the backreaction damps the decrease of the tension with 
respect to the damped fermion densities. There are also tiny regions near the axis, with areas inversely proportional 
to L e , and where T could grow with respect to one of the state parameters 1Z X or IZ^ (see Fig. [ll]). As previously, 
for reasonable values of /, the first discrete values of the parameters are out of these domains, and the tension always 
decreases with both fermion densities. Finally, the tension reaches negative values at large damped fermion densities 
(see Fig. U). 




FIG. 11. Variation domains of the tension, plotted in unit of M 2 , for spacelike currents with — 4ty < I < 0, in the 
(1Z X /M, TZ^/M) plane. These regions have the same geometrical properties as the line density energy ones in Fig. ^, but 
this time, the zones with vertical hatches are domains where T could grow with respect to 1Z X , whereas the horizontal ones 
correspond to growth with respect to TZ^,. For reasonable values of /, the discrete values of 1Z, represented by dots, are out of 
these regions, and the tension always decrease with both fermion densities. 
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FIG. 12. Sign of the tension for timelike currents with / > 0, and spacelike currents with —4n < I < 0, the curves have been 
plotted with L e = 10/M in aim at clearly separating the various regions. In the timelike case, the domain where the tension 
is positive is closed, the null tension curves intersecting the axes at 1Z = y2L e M 2 / 1 , whereas for spacelike currents, the null 
tension curves tends asymptotically to 72. = (27r/L e )\/2/(167r 2 — I 2 )- Recall that the regions where T < are unstable with 
respect to transverse string perturbations according to the macroscopic formalism. 



4- Relevant values of the parameters 



The previous derivation of the backreaction is built on the classical vortex background and it is acceptable only if 
the backreacted gauge fields do not perturb appreciably the Higgs and orthoradial gauge fields profiles (s ee F ig. ||). 
From Eq. (Jig), it will be the case only if Q l Qt and Q Z Q Z can be neglected compared to Q e Qg. From Eq. (133), and 

Q fl Qfl_.J*-U 



this condition reads 



(qc^) 2 ABZB 



QeQ 6 



K X K^ sin (6 X ~ f ) sin (9g - f) m 2 

TTTj 2 TTTj 2 



or as function of /(O^, O x ) and the damped fermion densities, 



AQEQ < 1, 



(149) 



j(e X) e*) 



TTT] 2 



< 1. 



(150) 



This condition is satisfied for damped fermion densities small compared to the string energy scale, or for tiny values of 
the function /(O^, B x ). Moreover, the backreacted gauge fields need to be small in order to no perturb signific antly 
the zero modes. From the equations of motion (|22|) and (|24|), this condition leads to AB, ££> <C TZ and from Eq. (133) 
to 



m b c Hx) 



E(A)Q(0) < 1. 



(151) 



On the other hand, the maximum value of / in Eq. (140) is clearly obtained when there are only particles or anti- 
particle s tra pped in the string (8 = o r 8 = 7r), and deriving the order of magnitude of the numerical integral 



in Eq. (141), using equations (134) and (|10q)-( 109 ) , one shows that the large values of I (as I > 47r) can only be 



obtained for model parameters which verify 



m b c i>(x) 
V 2 c<j> 



S(A)Q(0) > 1. 



(152) 



As a result, in order for the backreacted gauge fields not to modify the equations of motion of the fermions at first 
order, the function I has to be much smaller than 4-7T. If it is not the case, then the previous zero modes are no longer 
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(153) 



valid solutions and the relevant equations of motion, in the case of the fermions, now read, from Eq. ( p2| ) 

dr r 

^i + i 
dr r 

dr r 

^± + 1 
dr r 

where the angular dependence has not been written owing to Eq. ( |3l| ) and assuming mi = 7713 + n. The zero modes 
seem to acquire an effective mass proportional to AB or YjB. More precisely, they are no longer eigenstates of the 
7°7 3 operator since new spinor components appear [£1, £4 here, see Eq. (j3^) and Eq. (^2^ It is clearly a second 
order effect since the gauge coupling constant q can be removed in the previous equations (153) using Eq. (133) and 
assuming 
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with X the zero mode solution k = —uj for W fermions [see Eq. (|43)], and q 2 S£ the perturbation induced by the 
backreaction p5| . As a result, for strong backreaction, the semi-classical approach can no longer be used, since such 
second order effects appear as the semi-classical manifestations of the one loop quantum corrections, and thus, only a 
full quantum theory would be well defined. However, if there is only one kind of fermion trapped in the string, "J say, 
the zero modes are not affected by the backreaction since AB is only generated from the X current, and therefore 
vanishes [see Eq. (134)], so £ = X is always solution of the equations of motion (153), and identically for the X zero 



modes alone |34 . Note, that there is no contradiction with the usual index theorem since it is derived for Dirac 
operators, and thus without backreacted fields. This just shows that the modes propagating in the vortex with strong 
backreacted gauge fields are no longer well described by the usual zero modes. Physically, it might be the signature 
of a tunneling of the zero modes to another states. The massive modes which have not been considered here could be 
more relevant in such cases. 

On the other hand, the shape of the string might allow the fermion densities 1Z X and TL^ to reach the tiny regions 
where the energy decreases with one of them (see Fig. ||), by means of the zero mode vacuum quantum effects. The 
present toy model does not involve the effect of the radius of curvature R of the string, and it is reasonable that the 
contribution of the zero mode vacuum to the energy per unit length involves R through a redefinition of L e . If L e 
becomes smaller than L, the first discrete values of the fermion densities could be inside the hatched regions in Fig. ||, 
since the discrete values of the fermion densities only depend on the physical length of the string L. Note that it 
would therefore be necessary that the zero mode vacuum energy is negative, which is not the case without curvature 
in the simple framework of section IV B 1[ Such effects could be relevant for vorton stability, as, for a small radius of 
curvature, the string could become unstable to fermion condensation. 



Finally, the model can be used only at the tree order, and the conditions (150) and (151) are the validity criteria 
of the above derivations. 



VI. COMPARISON WITH THE SCALAR CASE 



Owing to the fermionic two-dimensional quantization along the string, the energy per unit length and the tension 
of a string carrying massless fermionic currents have been derived up to the first order in backreaction corrections. 
The state of the string is found to be well defined with four state parameters which are the densities of each fermion 
trapped in the string, and asymmetry angles between particles and anti-particles in each fermion family. It seems quite 
different than the bosonic charge carriers case, where the current magnitude is the only relevant state parameter [jlJJ , 
however, this is the result of the allowed purely classical approach where the superposition of many quantum states 
can be view as only one classical state owing to the bosonic nature of the charge carriers. As a result, there is a 
degeneracy between the number of bosons trapped in the string and the charge current. The quantization introduced 
to deal with fermions naturally leads to separate the charge current from the particle current through the existence of 
anti-particle exitations. Moreover, the magnitude of the current can only modify the equation of state at non-zeroth 
order because the chiral nature of fermions trapped in the string requires simultaneous exitations between the two 
families to lead to non- lightlike charge currents. 
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Nevertheless, some global comparisons can be made with the scalar case. First, for reasonable values of I <C 47r, the 
energy per unit length grows with the fermionic densities, whereas the tension decreases with them. However, note 
the relevant parameter for the change in behaviors of the tension and line density energy is the function I(Q X ,Q^) 
instead of the current magnitude in the scalar case. As it was said, / quantifies, through the asymmetry between 
the number of particles and anti-particles trapped in the string, the efficiency of the charge current per particle to 
be timelike or spacelike. The more positive is /, the more timelike the fermionic charge current per particle will be, 
and conversely the more negative I is, the more spacelike it will be. Once again, this difference with the scalar case 
appears as a result of the degeneracy breaking between particle current and charge current due to the fermionic nature 
of the charge carriers. 

The stability of the string with respect to transverse perturbations is given from the macroscopic formalism by the 
sign of the tension (for line density energy positive) |l9| ] , and we find that instabilities always occur for densities roughly 
close to M/y/An + 1, in finite domain for timelike current, in infinite one for spacelike currents with — 4tt < I < 0. 
Another new results are obtained from the multi-dimensional properties of the equation of state, in particular the 
problem of stability with respect to longitudinal perturbations differs from the scalar barotropic case where the 
longitudinal perturbations propagation speed is given by c\ = —dT/dU [ [Hill , and therefore its two-dimensional form 
has to be derived to conclude on these kinds of instabilities. Nevertheless, by analogy with the scalar case, since, in the 
non-perturbed case and in the infinite string limit, the equation of state verifies U + T = 2M 2 |^2| , the longitudinal 
perturbation propagation speed might be close to the speed of the light, even with small backreaction, and therefore, 
only transverse stability would be relevant in macroscopic string stability with massless fermionic currents. 

VII. CONCLUSION 

The energy per unit length and the tension of a cosmic string carrying fermionic massless currents were derived 
in the frame of the Witten model in the neutral limit. Contrary to bosonic charge carriers, the two-dimensional 
quantization required to deal with fermions, leads to more than one state parameter in order to yield a well-defined 
equation of state. They can be chosen, at zeroth order, as fermion densities trapped in the string regardless of charge 
conjugation. The minimal backreaction correction appears through the fermionic charge current magnitude which 
involves the asymmetry angles between the number of particles and anti-particles trapped in the string, and which 
might be identified with the baryonic number of the plasma in which the string was formed during the phase transition. 
As a result, it is shown that fermionic charge currents can be lightlike, spacelike as well as timelike. Moreover the line 
energy density and the tension evolve globally as in the bosonic charge carriers case, but it was found that the tension 
can take negative values in extreme regions where the fermion densities are close to the string mass, and where the 
string is therefore unstable with respect to transverse perturbations according to the macroscopic formalism. 

The present model has been built on the generic existence of fermionic zero modes in the string and follows only 
a semi-classical approach. It is no longer valid for higher corrections in the backreaction when they modify notably 
the vortex background and seem to give effective mass to the previous zero modes. It may be conjectured, at this 
stage, that in a full quantum theory, the quantum loop corrections give mass to the zero modes for high currents and 
consequently might lead to their decay by the mean of massive states. Only chiral charge currents could be stable on 
cosmic string carrying large fermionic massless currents in such a case. Another possible effect, relevant for vortons 
stability, may be expected for loops with small radius of curvature, by means of the vacuum effects which could render 
the loop unstable to fermion condensation. 

It will be interesting to quantify such modifications on the equation of state in future works, as the effects of 
worldsheet curvature, and the modification of the density line energy and tension by the massive bound states. The 
field of validity of the model could therefore be extended to higher energy scales which would be more relevant for 
vortons and string formation. 
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